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Abstract 

Background: As for other non-model species, genetic analyses in quail will benefit greatly from a higher marker 
density, now attainable thanks to the evolution of sequencing and genotyping technologies. Our objective was to 
obtain the first genome wide panel of Japanese quail SNP (Single Nucleotide Polymorphism) and to use it for the 
fine mapping of a QTL for a fear-related behaviour, namely tonic immobility, previously localized on Coturnix 
japonica chromosome 1. To this aim, two reduced representations of the genome were analysed through 
high-throughput 454 sequencing: AFLP (Amplified Fragment Length Polymorphism) fragments as representatives 
of genomic DNA, and EST (Expressed Sequence Tag) as representatives of the transcriptome. 

Results: The sequencing runs produced 399,189 and 1,106,762 sequence reads from cDNA and genomic 
fragments, respectively. They covered over 434 Mb of sequence in total and allowed us to detect 17,433 putative 
SNP. Among them, 384 were used to genotype two Advanced Intercross Lines (AIL) obtained from three quail lines 
differing for duration of tonic immobility. Despite the absence of genotyping for founder individuals in the analysis, 
the previously identified candidate region on chromosome 1 was refined and led to the identification of a 
candidate gene. 

Conclusions: These data confirm the efficiency of transcript and AFLP-sequencing for SNP discovery in a 
non-model species, and its application to the fine mapping of a complex trait. Our results reveal a significant 
association of duration of tonic immobility with a genomic region comprising the DMD (dystrophin) gene. Further 
characterization of this candidate gene is needed to decipher its putative role in tonic immobility in Coturnix. 
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Background 

The two avian Coturnix species - Japanese quail {Cotur- 
nix japonica) and Common quail {Coturnix coturnix) - 
belong to the order Galliformes and to the family 
Phasianidae, which also includes the chicken {Gallus 
gallus domesticus). The distributions of the two quail 
species overlap in the lake Baikal region, to and from 
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where wild quail migrate every year, toward Asia for C. 
japonica and toward Europe and Africa for C. coturnix 
[1]. However, hybridization is possible between the two 
species [2], and hybrid males have a distinct crow [3]. 

In Asia, the Japanese quail was domesticated for its 
singing abilities in the 15th century, and selected lines 
were developed in Japan for egg production at the begin- 
ning of the 20th century [4]. Japanese quail lines were 
then introduced successively in North America, Europe 
and the Near and Middle East in the first half of the 
20th century [5]. Its small size, short generation interval, 
good laying abilities and ease of reproduction interested 
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scientists looking for a bird easy to maintain and to ex- 
periment with in the laboratory, or willing to develop al- 
ternative ways of producing eggs and meat for human 
consumption [5]. However, for unknown reasons, the 
Common quail was not domesticated. Research on Japa- 
nese quail has led to the development and use of com- 
mercial quail strains for meat and egg production in 
Asia, Europe and other parts of the world [6] and to a 
better knowledge of its genetics [7] . 

The first modern genetic tools developed specifically 
for quail analyses were panels of markers [8,9] and gen- 
etic maps using AFLP (Amplified Fragment Length Poly- 
morphism) [10] and microsatellites [11]. Despite their 
limited extent, they have already been used extensively. 
At the population level, genome scans were carried out 
to detect QTL (Quantitive Trait Loci) for traits asso- 
ciated to fearfulness [12], growth, egg production and 
duration of tonic immobility [13] and egg laying curve 
[14], and markers were used to discriminate the two spe- 
cies and their hybrids [15,16]. At the genome level, genes 
carrying causal mutations were found for plumage col- 
ours, like TYRP1 for the roux plumage [17] and ASIP for 
the yellow [18] and recessive black [19] plumages. More 
recently, the existing quail maps were integrated and 
aligned to the assembled chicken sequence [20], thereby 
confirming the high degree of synteny conservation 
expected because of the high level of chromosome hom- 
ology observed between both species [21], and the rele- 
vance of quail models for trait analyses in poultry. 

Further developments in quail genomics are held back 
by the small number of available markers which delay 
the search for other causal mutations. It is especially 
critical for traits with no known homologous phenotypes 
in chicken, like the rusty plumage colour and the curly- 
shaped feathers [22], and for the fine mapping of QTL. 
Our objective was to develop a dedicated quail SNP 
(Single Nucleotide Polymorphism) panel, in order to fine 
map a QTL for duration of tonic immobility. Tonic im- 
mobility (TI) is an innate fear-related behaviour which 
occurs throughout the Animal Kingdom [23,24] and has 
been studied extensively in poultry (see [25]). Its dur- 
ation (dTI) is considered to be a measure of fearfulness, 
representing an antipredatory behaviour. In birds, it is 
characterized by a strong motor inhibition, elicited by a 
brief period of physical restraint. The animal first strug- 
gles and attempts to escape, then gets immobile even 
after termination of restraint, and this freezing-like pos- 
ture may last from a few seconds to several hours [24]. 
Duration of tonic immobility is a quantitative polygenic 
trait with low to moderate heritability in poultry [26]. It 
has been successfully selected upwards and downwards 
in experimental quail lines over a large number of gen- 
erations [27]. Despite its adaptive nature, fear may have 
deleterious effects on domestic animal welfare and health 



[25], and the relationship between duration of tonic im- 
mobility and productivity has been shown to vary from a 
positive to a negative one depending on the avian popula- 
tion and on the performance traits [28-30] . 

Deciphering the genetic architecture of tonic immobil- 
ity will be challenging, as for other complex traits [31]. 
We first aim to fine map a QTL associated with dTI pre- 
viously mapped to quail chromosome 1 by using two F2 
crosses between a Japanese quail line selected for its 
long tonic immobility (LTI) and two lines with shorter 
tonic immobility, STI (Short Tonic Immobility, selected 
for lower tonic immobility [12,13]) and DD (selected for 
early egg production [12,13]). To this end, two Advanced 
Intercross Lines (AIL) were built. Provided the density of 
information is sufficient, using advanced generations is 
expected to restrict the genome region involved in the 
variability of the trait [32], with the ultimate objective to 
identify a causative gene. We chose to obtain a high num- 
ber of SNP directly from our quail populations by high- 
throughput sequencing (Titanium 454 GS-FLX, Roche) of 
two main types of reduced representations of the genome: 
restriction digested fractions of genomic DNA, and EST 
(Expressed Sequence Tag), extending a first experiment 
[33] by increasing the number of sequences obtained and 
detecting reliable SNP. A subset of these markers was 
selected to analyse the AIL populations and to fine map 
the QTL. 

Results and discussion 

Sequence analyses and SNP discovery 

The sequences obtained in this study are available at 
NCBI (http://www.ncbi.nlm.nih.gov/, accession number 
SRP002189). 

Despite the rapid reduction of sequencing costs, identify- 
ing SNP by whole-genome sequencing of many individuals 
in a non-model species remains expensive, notably because 
the reliability of the identified SNP depends on the se- 
quence depth [34]. Therefore, another approach was 
chosen [35] and reduced representation libraries (RRL) 
were built: a genomic library with tagged samples (AFLP 
fragments) and a transcriptomic library (cDNA fragments). 

AFLP fragments were sequenced in one and a half run 
which produced 1,106,762 reads, totalizing about 280 
Mb of sequence. Reads were mapped to the chicken gen- 
ome (Figure 1) to determine if the chicken genome 
coverage was compatible with our aim to select quail 
markers evenly distributed across the genome, through 
the use of the chicken sequence assembly. The average 
number of reads per chromosome is consistent with the 
chicken chromosome lengths, with the exception of the 
Z chromosome, where more reads are observed than 
expected. This is probably due to the repetitive nature of 
a large part of the sequence of the Z chromosome [36]. 
In total, 72% of the sequences mapped to the chicken 
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Figure 1 Chicken genome coverage by quail AFLP fragments. A. Number of reads aligned to each chicken chromosome, after Burrows- 
Wheeler Alignment algorithm (bwa bwasw analysis). For each chicken chromosome (X axis), chromosome length in base pairs (bp, Chicken May 
2006 (WUGSC 2.1/galGal3) assembly) and number of mapped reads (count) is indicated (Y axis). B. Examples of chicken chromosome coverage by 
the quail reads for GGA01 (left) and GGA02 (right). X axis: position along the chromosome (physical length). Y axis: Number of reads aligned 

along the chromosome (counts, log scale), 
v J 



genome. This result is comparable with that obtained in 
turkey, using a different alignment algorithm [37], con- 
firming the possibility to obtain satisfactory results from 
alignment of heterologous sequences on a reference gen- 
ome assembly. This illustrates the fact that both species 
have a similar level of divergence (around 35 MYA) with 
the chicken [20,38]. Except for some regions with high 
depth due to repetitive elements, the AFLP coverage was 
uniform enough (Figure lb) to develop genome-wide 
SNP assays with confidence, given the proximity of the 
chicken and quail genomes [20]. 

From these data, 847,452 reads (76.6%) with identity 
tag and length higher than 120 bases were assembled 
into 42,860 contigs (N50 = 331 bp, consistent with the 
fraction of the AFLP fragments sequenced: between 200 
and 400 bp), leaving 95,812 singletons. In the absence of 
a reference quail genome sequence, the SNP detection 
was performed from the obtained contigs. This led to 
the identification of 10,545 SNP. 

Regarding the SNP type, a transition/transversion rate 
of 2.2/1 was observed. This value is similar to those 



reported for chicken (2.2/1), duck (2.3/1), and slightly 
higher than that for the great tit (1.7/1) [39,40]. Overall, 
our results confirm that reduced representation libraries 
from AFLP-derived fragments are successful tools for 
SNP detection. They also have been effectively used for 
phylogeographic studies [41-46]. 

In addition to the AFLP library, a second reduced rep- 
resentation of the genome was obtained from brain and 
embryo transcripts. The brain transcriptome may reflect 
a higher part of the genome than other tissues [47], and 
brain is probably the location of differential expression 
of behavioural genes between our quail lines. The 
embryo transcriptome was added to find SNP in gen- 
omic regions not covered by the brain transcriptome. 

As previously described [33], a total of 399,189 cDNA 
reads have been obtained, of which 374,939 filtered 
sequences were assembled into 31,010 EST contigs, with 
an average length of 733 bases (from 100 to 6393, N50 = 
786) and an average depth of 11 (from 2 to 3049), leav- 
ing 36,572 singletons. Among these 31,010 contigs, 
27,919 (90%) mapped to the chicken genome. 
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Out of the 67,582 contigs or singletons built from 
cDNA data, 54% had significant matches (evalue < 10" 5 ) 
on Swissprot, RefSeq_RNA or RefSeq_Protein databases. 
Sequences with more than 50% coverage with an identity 
level higher than 70% (n= 18,416) corresponded to 8,550 
unique referenced genes. 

As already established for several species including 
Brassica rapa, Drosophila melanogaster or Homo sapiens 
[48-50], the transition/transversion ratio is higher in 
coding sequences than in untranscribed regions, and the 
3.1:1 ratio in our EST data versus 2.2:1 from AFLP data 
confirmed this assessment. SNP were mainly discovered 
in non-coding region (Additional file 1: Figure SI), and 
30% of the observed SNP in the coding region were 
non-synonymous SNP, like in chicken [51]. 

A total of 17,433 SNP were detected: 10,545 came 
from the genomic fragments while 6,888 were identified 
from the cDNA sequences, 29 of which being common 
between the two analyses based on their mapping on the 
chicken genome. SNP with at least 25 bp of flanking se- 
quence each side have been submitted to dbSNP. This 
panel constitutes the first genome-wide collection of 
SNP markers for genetic studies in Coturnix. 

Tonic Immobility of AIL 

The distribution of dTI was similar for the two AIL 
populations (Additional file 2: Figure S2). Yet, the distri- 
bution for DD-LTI was shifted slightly toward higher 
values of dTI, and was therefore more affected by the 
censoring of dTI at 300 sec. This difference probably 
reflects the distinct origin and history of selection of the 
two lines with lower dTI, DD (not selected on dTI) and 
STI (selected for short dTI) which were used to start the 
two AIL populations in crosses with Line LTI (selected 
for long dTI). 

QTL analyses 

As previous experiments had shown the presence of a 
suggestive QTL on quail chromosome 1 [12,13], SNP 
were selected from sequences mapping to the homolo- 
gous chicken chromosome 1. A total of 1,098 SNP 
selected from AFLP and 1,567 observed in cDNA were 
obtained, and 384 SNP putatively covering the entire 
quail chromosome 1 were chosen for genotyping. 
Twenty-five of them, detected from quail cDNA puta- 
tively homologous to several chicken genes, appeared to 
map on different chromosomes when considering their 
best localization (Additional file 3: Table SI). Out of the 
384 markers, 352 (92%) gave reliable genotypes in our 
populations (AIL_DD: 326 SNP, AIL_STI: 336 SNP), and 
a total of 320 (91%) were polymorphic. After testing for 
non deviation from Hardy- Weinberg equilibrium, re- 
spectively 305 (AIL_DD) and 332 (AIL_STI) markers 
were used for further analyses. 



A framework map of quail chromosome 1 was built, 
comprising 60 markers with an average number of in- 
formative meioses of 505 (from 79 to 835), for a genetic 
length of 287.4 cM and with an average between-marker 
distance of 4.9 ± 3.9 cM (Figure 2). The length of our 
genetic map is comparable to that (274.8 cM) obtained 
previously for the same chromosome by Kayang and 
coworkers [20], but is shorter than the one (402.9 cM) 
obtained in a backcross study [52]. Several factors could 
explain this difference in length, like the structure of the 
populations and the number of individuals genotyped in 
each study, the type of markers used or the method 
employed to calculate the distance. 

Assuming quail and chicken chromosomes 1 have 
similar physical length, the relative coverage for our 
quail genetic map was 1.5 cM/Mb, which is much lower 
than the 2.4 cM/Mb ratio observed in chicken [53]. This 
result, probably partially related to the variability of re- 
combination rates between populations, is consistent 
with the smaller number of chiasmata reported in quail 
lampbrush chromosomes [54]. 

AIL populations are expected to offer a very powerful 
tool to fine map QTL, especially when inbred lines are 
used as founders [32], and the final AIL generation is ana- 
lyzed using classical QTL methods, taking advantage of 
the recombination events that have occurred over the 
course of multiple successive generations. We used this 
strategy with our lines to fine map a QTL previously 
detected in two distinct F2 crosses [12,13]. Two AIL 
populations were used: the AIL_DD families based upon 
F2 crosses of DD and LTI quail, and the AIL_STI popula- 
tion resulting from F2 crosses between STI and LTI birds. 

First, QTL detection was applied within AIL sire fam- 
ilies by interval mapping with QTLMap [55]. Unexpect- 
edly, no QTL for dTI could be detected in either of the 
final AIL generations, even after performing a bulk ana- 
lysis using the two AIL populations. We suspected this re- 
sult was the consequence of a lower detection power 
because the origin of alleles was not accounted for. In- 
deed, recent AIL studies included genotyping and pheno- 
typing of all individuals in the AIL line [56], or were set up 
from at least one inbred line [57], which made it possible 
to trace the origin of alleles back to the F0 lines. To over- 
come the problem, an association analysis was performed 
in our AIL, with the assumption that enough power would 
be gained by adding information from dams within fam- 
ilies of sires homozygous for the QTL, that were unin- 
formative in the QTL mapping analysis based on familial 
linkage. This strategy was successful in the AIL_DD popu- 
lation (Figure 3): 4 markers, in a genomic region ortholo- 
gous to position 121.54-121.78 Mb in chicken, showed a 
significant association with dTI (Bonferroni corrected for 
multiple testing p < 0.01): FSIEL0M02C8AXD_286, 
FSIEL0M02D00IT_258, FSIEL0M02C7RBX_180, and 



Fresard et al. BMC Genomics 2012, 13:551 
http://www.biomedcentral.eom/1 471 -21 64/1 3/551 



Page 5 of 1 1 



46,6- 
52,1 - 
54,2- 



287,4- 



F3SDQUK02CJYUN_291 

FSIEL0M02C111J_178 

FQU5R7M01 DVT66-0_330 

FQU5R7M01BZRAO-0_215 

FSIEL0M02D26QZ_39 

FQU5R7M02GF1 9M-0_272 

FQU5R7M01 AR2YZ-0_1 13 

FSIEL0M03GKREE_424 

FZVM8I102C1CRX_289 

FSIEL0M01BLHEM_251 

FSIEL0M02C2AQT_397 
FSIEL0M03GR4TM_261 
FSIEL0M02C6PJ3_261 

FSIEL0M01AZ4YN_263 

FZVM8I104I135K_332 

FQU5R7M02ITYZD-0_280 

FQU5R7M01A21EF-0_81 

FJXLD7H02H1NSI-0_231 

FSIEL0M02C45DP_108 

FQU5R7M01 CILD3-0_1 99 

FJXLD7H02ILSZS-0_82 

FSIEL0M03FL3NQ_267 

FQU5R7M01 CABVO-0_59 

FSIEL0M02C1CF5_155 

FQU5R7M02JKAQW-0_1 67 

FSIEL0M02DBPPH_27 



F3SDQUK02CDOVK_1 74 

FSIEL0M02C4EYCM38 

FSIEL0M02D2N49_373 

FSIEL0M02C9PVA_98 

FSIEL0M02EB774_354 

FSIEL0M02DF7NP_262 

F3SDQUK02B73C1_105 

FSIEL0M02D6BOC_162 

FSIEL0M02D06PR_377 

FQU5R7M01ADXMZ-0_197 

FQU5R7M02JJHRR-0_73 

FQU5R7M02I YSYE-2_1 42 

FQU5R7M01 D7NNZ-CM 86 

FSIEL0M01A1SZR_282 

FSIEL0M02C5W1M_237 

FSIEL0M02ERV8C_383 

FQU5R7M02F8IJS-0_1 59 

FSIEL0M01A4H2F_56 

FQU5R7M02FYIZK-2_1 37 

FSIEL0M02D5447_36 

FSIEL0M02DTOW8_461 

FSIEL0M03FTN66_287 

F3SDQUK02CIAEW_394 

FSIEL0M02D0JDA_268 

F3SDQUK02BU0JH_286 

FJXLD7H02HKLGL-0_1 62 

FQU5R7M02H140S-0_158 

FJXLD7H02JJLTX-0_66 

FSIEL0M02C9DBO_385 

F3SDQUK02B7324_256 

FQU5R7M01 B955I-CM 20 

FSIEL0M02C92LS_230 

FQU5R7M02H1ZJB-0_223 

F3SDQUK02BVBO3_470 



Figure 2 Genetic framework map of quail chromosome 1. 

Positions are given in cM. 



FSIEL0M02DMFKD_281 (Additional file 3: Table SI). 
Non-significant markers encompassing the region were 
respectively at 119.07 (F3SDQUK02B73C1_105) and 
123.28 Mb (FSIEL0M01A1W48_163). As the AIL ap- 
proach allows us to reasonably consider that linkage dis- 
equilibrium is limited, our investigation could be first 
restricted around this 4.21 Mb area. This result confirms 
the first mapping of a suggestive QTL in the original F2 
population, as the region detected here belongs to the 
confidence interval estimated by the two-LOD drop-off 
method [58] from the original LTI-DD QTL position 
[12,13,58]. The same region also contained a significant 
QTL in the original LTI-STI F2 population [12], but had 
not been previously shown to be associated with dTI in 
chicken [59]. 

The association study showed also that the most sig- 
nificant marker (FSIEL0M02C8AXD_286, orthologous 
position 121.54 Mb on the chicken genome) had a sub- 
stitution effect of 0.45 phenotypic standard deviations 
and that only 3 F6 sires were heterozygous at the QTL. 
By using these estimates, the power of detection of the 
QTL in the AIL_DD familial analysis was estimated to 
reach only 40%, which accounts for the fact that the 
QTL was not detected using this approach. 

The same analysis in the AIL_STI population did not 
reveal any significant association between dTI and a 
marker. This may be due to the lack of linkage disequi- 
librium between the markers and the causative mutation 
in the STI population: whereas the most significant mar- 
kers from the AIL_DD were fixed in the DD line (all F0 
parents homozygous for the same allele), the STI popu- 
lation showed allele sharing, indicating that a complete 
association between a putative fixed mutation and a SNP 
marker could not be observed in the AIL_STI. The facts 
that the AIL were not phenotyped under the same set- 
ting (housing, technician) and that the genetic origins of 
the F0 were a same base line population for LTI and 
STI, and a different unrelated one for DD [60] may have 
also contributed to the diverging results of the associ- 
ation analyses. 

Most of the genes recorded in the investigated region 
have no known association with tonic immobility 
(Additional file 4: Table S2), with the notable exception 
of the Dystrophin gene (Q9PTS7_CHICK or DMD). In 
the candidate region, one gene (IL1RAPL1, IL-1 -receptor 
accessory protein-like 1, 119.54-120.16 Mb) is involved 
in mental retardation or autism [61,62]; structural varia- 
tions including another gene (PTCHD1, patched domain 
containing 1, 121.89-121.94 Mb) have been shown to be 
involved in mental disorders like autism [63], but no 
clear connection could be made with tonic immobility. 
The 4 significant markers are included in 4 different 
genes (Additional file 4: Table S2): Q5ZLT2, or PDK3 
(pyruvate dehydrogenase kinase, isoenzyme 3), a kinase 
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involved in the regulation of glucose metabolism, has 
been showed to be associated with drug resistance and 
colon cancer [64]; EIF2S3 (eukaryotic translation initi- 
ation factor 2, subunit 3 gamma) is a gene coding for a 
subunit of eukaryotic translation initiation factor 2, 
involved in protein synthesis by playing a role in transla- 
tion start site selection (see [65]); KLHL15 (kelch-like 
protein 15) encodes for a protein belonging to the 
Kelch-like proteins family, implicated in embryogenesis 
and carcinogenesis mainly through cytoskeleton 
organization [66]; all these three genes do not seem to 
have a function directly linked with the duration of tonic 
immobility. The fourth one, SAT1 (spermidine/spermine 
Nl-acetyltransferase 1), catalyzing the acetylation of 
polyamines, has been associated with antineoplastic ac- 
tivity of several antitumor [67] and is essential for the 
polyamine homeostasis (see [68]). SAT1 has been shown 
in transgenic mice to be involved in several behavioural 
traits as activity and spatial learning in females, but not 
in defensive reactivity [69]. It may constitute a correct 
candidate gene. But the DMD gene, localized between 
118.07 and 119.07 Mb, seems to be a very interesting 
one even if only the 3 end of the gene is included in the 
candidate region. Sekiguchi et al. [70] showed that DMD 
is involved in defensive behaviors in mouse, including a 
strong freezing response, reminiscent of the strong 
motor inhibition exhibited by quail in the tonic immo- 
bility response. This gene, well known for its implication 
in Duchenne muscular dystrophy in Human through 
mutations leading to loss-of-function of the gene products 



in striated muscles [71], is also expressed in the brain, 
where it might account for the mental deficiency reported 
in Duchenne muscular dystrophy patients (see [72] for a 
review). However, further work is needed to investigate 
the precise involvement of brain dystrophin in the control 
of emotions and cognition. 

Conclusions 

Using Japanese quail as an example, this work confirms 
the efficiency of AFLP fragments as a reduced represen- 
tation of the genome to develop SNP markers for spe- 
cies for which no whole-genome draft sequence is 
available. Transcript sequencing brings a high number 
of quail EST that can be used to help annotate a future 
quail genome assembly. This study provides also a high 
number of Japanese quail SNP, for genetic analyses in 
Coturnix. It underlines the importance of modern tech- 
nologies as Next Generation Sequencing in the renewed 
ability to perform association analyses, even in non- 
model species. 

The association study carried out with a SNP subpanel 
for chromosome 1 led to a candidate gene, DMD, as 
being partly responsible for the variability of tonic im- 
mobility in Japanese quail. It should be the focus of fur- 
ther investigations to help unravel this complex trait. 

Methods 

Experimental animals 

Tonic immobility was the selection criterion for the di- 
vergent selection of two quail lines: the LTI and STI 
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lines [27]. These two lines and a third one, line DD 
unselected on tonic immobility and of a different origin 
[73], were used for the study. Two AIL populations were 
produced: the "AIL_DD" population was obtained by 
crossing 6 DD male to 6 LTI females and 6 LTI males to 
6 DD females to produce the Fl generation. The off- 
spring were then intermated within each generation 
until the F7 generation (603 offspring, Table 1). The 
n AIL_STI" population was founded by crossing LTI and 
STI quails and carried out until the F6 generation (679 
offspring, Table 1). Between 9 and 11 days of age, dur- 
ation of tonic immobility was measured on each quail of 
the final AIL generation. TI duration is the time during 
which a bird remains immobile after the fear reaction has 
been induced by keeping it on its back in a U-shaped cra- 
dle for 10 seconds (Additional file 5: Figure S3). dTI was 
set to 0 second if tonic immobility was not induced after 
five attempts. If the quail was still in a tonic immobility 
state after 5 min, the observation was censored and dTI 
set to 300 seconds. Statistical analyses were run on the 
transformed variable LOG(dTI+l). 

454 sequencing 

Transcripts and AFLP fragments libraries' preparation, 
and 454 sequencing, have already been described [33]. 
Three quail lines were used: LTI, STI, DD [27,73], and 
total RNA was extracted from about 500 mg of adult 
quail brains (4 samples from each line) and total 
embryos (Embryonic day 8, 2 samples from each line) 
for library preparation. Blood samples from two indivi- 
duals from each line were used for the AFLP library 
preparation. 

Sequence analyses 

Trimming of low-quality sequence and 454 adaptors were 
performed using the Newbler software (Roche/454). 

AFLP analyses 

Identification of tags within each sequence was per- 
formed by using crossmatch (from the PhredPhrap suite, 
http://www.phrap.org/phredphrapconsed.html, Phil Green) 
with tagged AFLP primers as query sequences. The 
crossmatch results were used to split the initial sequences 
when primers were observed within the fragment se- 
quence (and not only at its ends). We then selected the 
sequences long enough for the clustering step (>120 bp). 
Sequences with 1 or 2 tags were used to generate se- 
quence clusters, based on more than 95% identity (over a 
length of 80% of the shorter sequence) between cluster 
fragments, as processed by Megablast [74]. Next, contigs 
were assembled from the clusters with the CAP3 sequence 
assembly program [75]. 

SNP detection was performed with an in-house script: 
we selected SNP with a minimum depth of 5 sequences, 



a minus allele count of 2 and having complete similarity 
between the sequences in the ten base pairs located on 
the left and right sides of the SNP. 

Alignment of the reads to the chicken genomic se- 
quence (Chicken May 2006 (WUGSC 2.1/galGal3) assem- 
bly) was performed with the Burrows-Wheeler Alignment 
algorithm bwa, version 0.6.1-rl04, command bwasw, at 
default parameters (mismatch penalty of 3, maximum one 
gap opened, gap open penalty of 11, gap extension penalty 
of 4. . .). This algorithm, more permissive about alignment 
gaps, has been developed for long queries but was used 
here to deal with heterologous alignment [76]. Genome 
coverage was calculated from the SAMtools v0.1.18 
(mpileup tool) utilities [77] (samtools mpileup -A -B) 
and visualised through a custom script. 

cDNA analyses 

Artefactual sequence duplicates due to the 454 technol- 
ogy (98% identity, same start, same strand, same 454 
run) were removed [78]. 

Because of the large number of sequences, we used a 
two-step process to assemble these sequences into contigs. 
The first step built clusters of sequences sharing at least 
75 bp with an identity rate of 96% using MegaBlast [74]. 
The second step constructed coherent contigs from the 
previous clusters using CAP3 [75], at the recommended 
stringency of 40 bp overlap with 90% sequence identity. 

As no annotated quail reference genome is available, the 
contigs or singletons were annotated by BLAST [79] on the 
chicken genome using several databases (SwissProt, 
Refseq_RNA, RefSeq_Prot). Sequences with more than 50% 
coverage and with an identity level higher than 70% were 
considered as homologous to a reference gene, and their 
GO annotation were retrieved from the databases. All data 
were loaded in a locally-adapted Ensembl database. 

The SNP detection was performed with the same script 
and parameters as for the AFLP fragments. Functional 
annotation of the SNP was performed with custom Perl 



Table 1 Number of animals in each generation of the AIL 
pedigrees 



Generation 




AIL_DD 






AIL_STI 




Offspring 


Sires 


Dams 


Offspring 


Sires 


Dams 


F0 


24 


/ 


/ 


72 


/ 


/ 


Fl 


40 


12 


12 


72 


36 


36 


F2 


433 


10 


26 


72 


36 


36 


F3 


96 


14 


38 


72 


35 


35 


F4 


98 


48 


48 


72 


36 


36 


F5 


98 


48 


48 


71 


32 


32 


F6 


52 


35 


35 


679 


10 


30 


F7 


603 


13 


39 


/ 


/ 


/ 
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scripts, using the Ensembl APIs (Application Programme 
Interface) [80]. 

Marker genotyping 

In order to select SNP from contigs uniquely mapped to 
GGA01, a Megablast analysis [74] was performed to 
map the AFLP contigs on the chicken reference genome 
(Chicken May 2006 (WUGSC 2.1/galGal3) assembly). 
Similarly, cDNA SNP were selected from genes localized 
on chicken chromosome 1. Then, by using their ortholo- 
gous positions on the chicken genome, 384 SNP were 
chosen as putatively covering the entire quail chromo- 
some 1 (Additional file 3: Table SI). They were analysed 
on 1399 animals (AIL F6 or F7, their parents, and 26 FO 
individuals, corresponding to the 24 AIL_DD founders 
and 2 LTI individuals) using the Illumina GoldenGate 
assay. Genotype calling was performed with GenomeStu- 
dio V2010.1 (Illumina). Genotypes were filtered for call 
rate (> 0.85), call frequency (> 0.50) and cluster disper- 
sion. Markers with an unexpected deviation from 
Hardy- Weinberg equilibrium (p < 10~ 6 ) were removed 
with PLINK (option - hardy, http://pngu.mgh.harvard. 
edu/purcell/plink/ ) [81]. 

QTL analyses 
Interval mapping 

A genetic map of quail chromosome 1 was built by using 
the CRI-MAP program [82] version 2.503 (revised ver- 
sion of Phil Green's CRI-MAP v2.4, modified by Jill 
Maddox and Ian Evans) with the AIL_STI population, 
which had a higher number of reliable markers. The 
"twopoint" option allowed to group markers in a unique 
linkage group (twopoint LOD score >= 3). The "build" 
option was used to order markers within the linkage 
group, while the "flips" option enabled to confirm the 
order of the markers. We build a framework map with a 
LOD threshold of 3, using a stepwise locus-adding strat- 
egy, starting from the markers with the higher number 
of informative meioses. 

QTL interval mapping was performed with the 
QTLMap software using likelihood computation at every 
centimorgan, ignoring the maternal chromosomes [55]. 
The linearized likelihood under the alternative hypoth- 
esis was computed as: 

L x =Z i Z k f\y ik /hs h M i \, 

with f\y ik /hs h M^N \Z p b p a p + yi t + Pfca?] a t ) , 

with i the sire family in the design; ik the progeny within 
sire family i; j\yiklhs u M^ the likelihood of individual ik 
phenotype computed conditionally to hs t the sire i 
phase and M t the marker information in the sire family 
i, considered as normally distributed with mean com- 
prising some nuisance fixed effects a p) the sire family 



mean of the phenotypes {i u and the QTL substitution 
effect in sire family i a*, and standard deviation of the 
sire family i o t (heteroscedastic model). P\ k is the trans- 
mission probability of one of the sire chromosome to its 
progeny ik at position x. The technician having recorded 
the phenotypes, the sex and the hatch were set as nuis- 
ance fixed effects in the model. Chromosome-wide sig- 
nificance thresholds at 5% were calculated through 1,000 
phenotype permutations. Power of detection was tested 
on AIL_DD modulating its two principal components, 
QTL frequency (ie proportion of heterozygote sires in 
the parental F5 generation) and QTL effect, by simulat- 
ing QTL using QTLMap. 

Association analyses 

Marker-trait association was tested by using a linear 
mixed model: 

y = Xfi + Z 1 a + Z 2 b + e 

where 

P : fixed effects (intercept, hatch, sex, technician, geno- 
type at marker) 
a : random sire effect on trait mean a~N(0,Io*a) 
b : random dam effect within sire on trait mean 

b~N(0jag) 

e : residual e~N(0,Io*^) 
where X, Z 2 and Z 2 are design matrices and / denotes 
the identity matrix. 

The model was implemented using the R « nlme » 
package {Linear and Nonlinear Mixed Effects Models) 
[83]. The Bonferroni correction was applied to correct 
for multiple testing. 

Candidate gene identification 

A search for candidate genes was performed through 
AnnotQTL [84] (http://annotqtl.genouest.org). This soft- 
ware gathers the functional annotation of genes mapped 
to a specific chromosomal region, from several public 
databases (Gene Ontology, Mammalian Phenotype, 
HGNC and Pubmed), in order to help find the best can- 
didate genes. 

Additional files 



Additional file 1: Figure SI. Distribution of the functional annotation 
of the detected SNP (log scale). Annotations are from the Ensembl APIs 
(Application Programme Interface): the functional consequence of each 
SNP in each transcript has been predicted using the Variant Effect 
Predictor (VEP). All the non intergenic consequences are represented. 

Additional file 2: Figure S2. Distribution of the trait in both AIL 
populations. Tl : Tonic Immobility. 

Additional file 3: Table SI. List of the 384 SNP used in this study. 
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Additional file 4: Table S2. List of the genes localized in the candidate 
interval. Source : Biomart, Ensembl Genes 61 (Sanger UK), Gallus Gallus 
Genes (WASHUCS52). 

Additional file 5: Figure S3. Quail in a tonic immobility state. 
Immobilization is induced by keeping the animal on its back in a 
U-shaped cradle for 10 seconds. 
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